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1. As is well known, the renormalization group (RG) is a tool providing 
deep insight into the structure of a quantum field theory. It is certainly worth 
applying this tool to quantum gravity. This work is devoted to the develop- 
ment and the application of the real space renormalization group technique 
in the context of euclidean quantum gravity. It is a direct continuation of 
the work 

We choose the lattice gravity framework. More precisely we adopt the 
particularly promising dynamical triangulation approach . The remarkable 
results obtained within a class of exactly solvable models in two dimensions 
strongly suggest that the dynamical triangulation recipe is the correct way 
of discretizing gravity (at least for fixed topology). 

In conventional statistical mechanics a real space renormalization group 
transformation has two facets: 

(a) geometry - cells of the body are "blocked" together. 

(b) matter fields - "block" fields are defined in terms of the original fields. 
On a regular lattice it is trivial to perform the step (a) in such a manner 

that the resulting lattice is identical, modulo rescaling, to the original one. 
Since the values of critical couplings depend on the lattice type this self- 
similarity feature of the transformation is important. On a random lattice 
an appropriate definition of (a) requires some thought. In this work we con- 
sider pure geometry, without matter fields, and consequently the geometrical 
aspect of the renormalization group. 

In ref. [Q a method of "blocking" triangulations that exploits the self- 
similarity feature of random manifolds has been proposed. Without repeating 
in detail the arguments presented in U let us briefly sketch the main idea . 
The intuitive arguments given below will be replaced progressively by more 
precise ones later on. We do not wish to give an impression of complexity 
from the outset. 

2. In 2d one can show J| that an infinite randomly triangulated manifold 
is a self-similar tree obtained by gluing together sub-structures called baby 
universes (BUs), defined as subuniverses separated from the remaining part 
of the universe by a narrow neck. We conjecture that a similar picture holds 
in 4d, at least in the neighborhood of the phase transition point. The results 
presented in this paper strongly support this conjecture. 

3 As we have learned from J. Ambj0rn during the LATTICE '94 Conference, our ideas 
partly overlap with those discussed earlier for 2d in ref. H . 
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We have proposed U to define the step (a) as the operation of cutting 
the last generation of BUs, i.e. those which have no further BUs growing on 
them. When the tree is finite, and as one continues cutting the successive 
generations of BUs, it gets smaller (in lattice units) and looks less and less 
branched. In the spirit of the renormalization group this is interpreted as 
reflecting the reduction of the resolving power: although, in physical units, 
the manifold remains of the same size, the fine structure is not observed. 

The next problem consists in establishing a connection between lattice 
and physical observables. For example, in lattice QCD one identifies the 
inverse of the correlation length £ with the mass m of a hadron 

ma = £ _1 (1) 

Here a denotes the lattice spacing and £ = £(«) depends on a dimensionless 
coupling constant k. The correlation length £ can be calculated for a given 
k. The mass m is a physical observable, independent of the regularization 
procedure. Hence, eq. (|1|) yields the variation of the lattice spacing a, mea- 
sured in physical units, with the coupling k. The so-called /^-function can by 
defined as 

where the derivative is taken along the RG flow. In QCD, the lattice cor- 
relation length diverges exponentially : £ oc exp(l/2/? K 2 ) with 
Po > 0. Eq. (JT|) implies that (3(k) = — (3qh?. The fixed point is at k = 
and is ultraviolet (UV) stable. The lattice spacing tends to zero, in physical 
units set by eq. ([!]), as one approaches the fixed point. 

We recall these very well known things, in a sense, to set the standards. 
We would like to have similar arguments in quantum gravity. The phase 
transition observed in 4d P, |j seems to be continuous and therefore one has 
a candidate for a fixed point. The nature of this fixed point (IR or UV) is to 
be found. The soundness of the proposed renormalization group has to be 
further examined. Finally, one has to select appropriate observables in order 
to fix the scales. We shall report about some progress achieved in studying 
these problems. But before proceeding we have to define the theory we are 
working with. 

4. We take an euclidean version of the Einstein-Hilbert action which for 
a simplicial manifold reads 

S = -k 2 N 2 + k 4 N 4 , (3) 
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where N 2 and N4 denote the total number of triangles and 4-simplexes, re- 
spectively. The theory is defined by the partition function 

Z(k 2 , Ka) = E Z N2Ni e~ s (4) 

N 2 ,N 4 

where 

Zn 2 n 4 = E W(T) (5) 

T(N 2 ,N 4 ) 

The sum in eq. (§) runs over fixed topology 4d simplicial manifolds T(N 2 , N&) 
with N 2 triangles and N4 4-simplexes. W(T) is the symmetry factor taking 
care of the equivalent relabelings of the manifold. 

We study a fixed volume N4 = N canonical ensemble for spherical mani- 
folds with the partition function : 

Z(k 2 ,N 4 ) =J2Z N2N4 e- s (6) 

N 2 

In practice, in order to have an ergodic algorithm, one is forced to simulate 
a multicanonical ensemble (see e.g. [0). One lets the volume fluctuate in an 
external potential U(N 4 ) with parameters adjusted in such a way that the 
entry N4 = N, at which measurements are done, is highly probable in the 
resulting N4 distribution. In our simulations we have used U = |(iV — N4) 2 
with 5 = 0.001. 

The observables have to be geometrical objects. One obvious choice is 
the total physical volume V = Na 4 . Once V has been fixed, the length of 
the lattice step in physical units is also fixed for a given lattice: 

a=(V/N) 1 * (7) 

However, the continuum limit of the theory is not yet defined. A second 
observable is needed for that. A priori, this second observable could be 
the local curvature R. It is, however, not quite obvious how the physical 
curvature should actually be defined. In numerical simulations of dynamically 
triangulated 4d manifolds, the bare Regge curvature residing on the hinges 
of the lattice is found to be nonvanishing at the transition point ^, g], 
which leads to basic problems in defining its dimensionful counterpart in the 
continuum limit. Various interpretations have been given to this finding. We 
are tempted to adopt here the point of view formulated by the authors of ref. 
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M, who argue in substance that the bare curvature is a lattice artifact and 
is not the curvature that would measure an observer using a physical stick, 
independent of the cut-off. 

Following the suggestion put forward in ||, we choose as our second 
observable a length scale characterizing the large-distance behavior of the 
so-called puncture-puncture correlations function, viz. the average geodesic 
distance between two randomly selected simplexes of the manifold: 

(r) = (N- 2 Y,ra b )N , (8) 

ab 

where r a b is the geodesic distance between 4-simplexes a and b and (...) n de- 
notes the average over manifolds described by the canonical partition function 

In practice, it is difficult to identify BUs with necks larger than the min- 
imum one. Following the terminology introduced in ref. [|[] we call minBUs 
the baby universes with the smallest possible neck. We have sugested in [|TJ 
that a meaningful RG transformation is defined by the operation of cutting 
last generation minBUs. Using the pure 2d quantum gravity as a toy model, 
we have checked that, on the average, the area of surfaces obtained after 
cutting last generation minBUs is, for large surfaces, a constant fraction of 
the area of the original surface. It is important to check that a similar result 
holds for (r), if (r) is to be used in a RG analysis. 

5. We have been working with surfaces of area ranging between N = 4096 
and 65536 triangles. The minimum number of triangles in a minBU has been 
set to 15. The number of heating sweeps between measures has been set to 
twice the autocorrelation length determined, for the observable (r), in ref. 
1(J Q. The results are shown in Table 1. The subscripts "in" and "out" refer 



to the original and the cut surfaces, respectively. Actually, the size of the 
original surface has been kept fixed : Ni n = N. 

It appears from the results presented in the Table 1 that the ratios Xn = 
{N out) /Nin and X r = (r out ) / (r in ) converge towards limiting values when N — > 
oo. Furthermore, the critical exponent v defined by the relation A r = \ V N 
approaches the expected value \. Let us recall that the problem of scaling in 

4 Strictly speaking, they have done this for a model with 2d gravity coupled to a gaussian 
field; it does not seem that autocorrelation lengths are longer for pure gravity when one 
measures such a purely geometrical object as the average geodesic distance. 
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pure 2d gravity has been solved rigorously |TT[| and it is known that (r) oc 
as N — > oo. This could have been checked numerically with better precision 
using a more direct method (see e.g. ref. Our point is that (r) scales 

correctly under the RG transformation we have proposed. Notice, that self- 
similarity and scaling only hold in the statistical sense: it is the grand- 
canonical ensemble that is mapped into itself by our transformation. 
In 2d gravity one can write for small enough a 

(r) = (p/ar (y = \) (9) 

where p is an observable with the dimension of a length, which indirectly 
determines the physical curvature of the surface. Eq. (|9|) is an analogue of 
eq. (0) and p could be used to set the scale of the theory. However, once the 
physical area V = Na 2 has been used for that purpose, the choice of p is no 
longer free. 

Playing with the toy model we have convinced ourselves that (r) is an 
observable suitable for our purposes. We also got better acquainted with the 
proposed RG transformation. Thus, when one cuts the successive generations 
of minBUs, the corresponding scaling factors A^, , ... , ... are not at 
all equal and actually rapidly tend to unity. This is easily understood from 
the results of ref. jp: before the outer minBUs have been cut, the minBUs 
of the next layer are obviously more volumineous than the outer ones (since 
the latter are just outgrows off the former). In order to keep A at constant, 
one should enlarge the necks of BUs to be cut at each iteration of the RG 
transformation. This is not what we are doing, since we restrict ourselves, 
for purely technical reasons, to minBUs. 

The results shown in Table 1 have been obtained iterating the transfor- 
mation until no minBU is left. When only one iteration is performed, the 
value of Aa? is larger (about 0.86 , with the same definition of a minBU) and 
A r is very close to unity. Consequently, the error in the determination of v is 
much larger, for the same statistics, and the results appear less convincing. 

6. Encouraged by the results obtained in 2d we go over to 4d. Now, (r) 
is a function of two variables, viz. N and k>2, to be denoted k from here on 
for the sake of simplicity of writing. 

Suppose we start with iV = N in and k = K in and calculate the corre- 
sponding average geodesic distance (r) = (r in ). After one iteration of RG we 
find (N out ) and (r out ). Denote by K out the value of k characterizing the new 
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ensemble. Define 



From (0) we get 



Sr = (r in ) - (r out ) (10) 
5N = N in -(N out ) (11) 
5k = Ki n K ou t (12) 



<51n(l/a) = ^S\nN (13) 
The associated shift of k can be found from the equation 

5r = r N 5\n N + r K 5n (14) 

where and r K are the partial derivatives of (r) with respect to IniV and 
k, which can also be estimated from numerical simulations. Once 8k has 
been calculated, the /3-function is readily found from (^) and flT3|) . Writing 
eq. (15) we have actually conjectured that k is the only relevant coupling 
at long distances. This is presumably reasonable when k is close enough to 
its critical value. One can check the conjecture by calculating the /3-function 
using some other observable than (r), to fix the scale. If the conjecture is 
true the result should not depend appreciably on the choice of the observable. 

We have applied our RG transformation to 4d spherical simplicial com- 
plexes with 8000 simplexes, working with a series of values of k ranging 
between 0.9 and 1.45 . The average geodesic distance has also been calcu- 
lated, for the same values of k at N = 6000, in order to estimate r^- The five 
standard local moves have been used for updating. Typically, we have been 
carrying ~ 10 5 sweeps of the lattice, for each k. Close to the critical point 
the number of sweeps has been significantly larger. Measures were separated 
by a number of heating sweeps, equal to twice the autocorrelation length, 
estimated in advance. The errors were found using the conventional binning 
method. The minimum volume of a minBU has been set to 20 simplexes. 
The RG transformation has been iterated only once. Above k = 1.15 cutting 
more than one generation of minBUs yields too small (N out ) for the differen- 
tial formula ([14]) to be trusted The data are collected in Table 2 and the 



3 At n = 1.3 , iterating RG the maximum number of times one finds (N ou t) = 859±300: 
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/^-function is shown in Fig. 1. A linear fit to the data yields /3(k) = (3q{h* — k) 
with f3 = 6.5(1.2) and k* = 1.2(2). 

Although the errors are large, especially as one enters the branched poly- 
mer phase, where the local moves are not very efficient, it is rather evident 
from Fig. 1 that the critical point is an ultraviolet stable one^|. Notice, that a 
qualitatively identical result has been found analytically in 2 + e dimensions 
[ |I"3l , IHH , using perturbation theory together with the e expansion. Let us now 
briefly outline some general consequences of this result. 

7. Close enough to the fixed point k = k* one can write 

lm - r M^- K) (/%>(.) (15) 

Solving this differential equation one finds 

a = a | n* - k | 1//3 ° (A)>0) (16) 

where ao is an integration constant, analogous to A -1 in QCD, and which 
should be given a value, in physical units, in order to define the theory. 
Comparing with (0) one gets 

V/4 = t (17) 

with t defined by 

t = N | k* - k | 4/A) (18) 

The quantity on the left-hand side of (|TT| ) is a constant, determined by 
the choice of V and ao- The trajectories of the RG flow are the curves t = 
const. The continuum limit is the double-scaling limit: iV — > oo, k — > n* 
and t = const. In analogy to 2d one can write (r) = (p/a) 4u , but now p and 
v may depend on the value given to t. 

8. Going farther ahead is beyond the scope of this paper. A few remarks 
are in order at this point : 

(i) This study should and will be extended to larger lattices. In order to do 
that it will be necessary to use a more efficient updating algorithm. Indeed, 

the manifolds resemble branched polymers and the situation differs from that described 
in sect. 5. 

6 P. Bialas has informed us that using the data of ref || and our definition of the 
/3-function one can estimate the latter in 3 points. The estimates agree with our result. 
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the local algorithm becomes inefficient as one enters the "cold" phase of 
the theory, where the manifold develops a branched polymer structure. The 
"BU surgery" algorithm proposed in |TU|, [T5| should be appropriate there. 
Performing more intensive MC studies one should also be able to estimate 
the corrections to the formula (|l~4"D coming from irrelevant directions at the 
fixed point, neglected here. 

(ii) Remember, that v~ x has the significance of the intrinsic Haussdorf 



dimension dn of the manifold ||15|| . It would be interesting to look for the t- 
dependence of du and determine the class of theories satisfying the constraint 
dn = 4. Our experience with 2d indicates that this may require going to very 
large systems. 

(iii) It would be very interesting to extend the analysis to models where 
matter fields are coupled to gravity and to achieve a better understanding of 
the relation between the results generated by lattice models and the contin- 
uum calculations in 2 + e dimensions. With a more efficient algorithm and 
a better statistics it will hopefully be possible to examine the dependence of 
the shape of the /5-function on the number of matter fields. Calculations in 
the continuum [13|, |14]] suggest something analogous to the c = 1 barrier of 
2d models: the number of matter fields cannot be arbitrary for the theory 
to exist. It would be important to obtain such a result in lattice gravity, 
without any recourse to perturbation theory. It might be that this has some 
relevance for the problem of the number of generations. 

In summary, we have applied the previously proposed real space renor- 
malization group procedure to 2d and 4d models of pure euclidean quantum 
gravity. In 4d we have defined a /^-function and we have found numerically 
some evidence for the existence of an ultraviolet stable fixed point of the 
theory. We have concluded mentioning a few problems for the future. 
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Figure caption 

Fig. 1 - The /3-function versus the coupling k in 4d, calculated for spherical 
manifolds with 8000 simplexes. 
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Table 1: Summary of data in 2d. 
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Table 2: Summary of data in 4d. 
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